The spin-Peierls transition beyond the adiabatic approximation 
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We develop a theory of the spin-Peierls transition taking into account the three dimensional 
character of the phonon field. Our approach does not rely on the adiabatic or mean field treatment 
for the phonons. It is instead based in the exact integration of the phonon field, the exact long 
wavelength solution of the one chain spin problem, and then a mean field approximation for the 
interchain interaction. We show that the spin gap and the critical temperature are strongly reduced 
due to the finite frequency effects of the phonon coupling transverse to the magnetic chains. We 
claim that our results should be applicable to the inorganic spin-Peierls compound CuGeOs. We 
show that the long standing discussion on absence of a soft mode in this compound can be naturally 
resolved within our theory. 

PACS numbers: 75.10.Dg, 75.10.Jm, 63.70.+h 

I. INTRODUCTION 

Quasi-one-dimensional electronic and magnetic materials are fascinating objects which have been intensively studied 
during the last decades. They show a variety of different phases as the magnetically dimerized, unconventional metal 
and triplet superconductor phases. It is possible to go from one phase to another by changing the experimentally 
accessible macroscopic parameters: pressure, temperature and the applied magnetic field. The different phases arise 
from the competition between the Coulomb repulsion between electrons and their coupling with the lattice degrees of 
freedom. Whereas the essential physics relics on the one-dimensionality of the electronic dynamics, the stabilization of 
the different phases depends on the three dimensional character of the system and thus on the interchain interaction. 
Transitions between the different phases are, in a way, a dimensional crossover from ID to the three dimensionality 
of any real material. 

Spin-Peierls systems are probably the simplest example of quasi-one-dimensional materials. They show a high 
temperature phase which can be well modelled by a one dimensional antiferromagnetic Heisenberg chain and a low 
temperature phase with lattice dimerization and spin singlet ground state. The transition between them arises from 
the competition between the gain in antiferromagnetic energy and the loss in elastic energy. It is indeed possible 
due to the three dimensional character of the phonon field. It is therefore very attractive to carefully understand the 
precise character of this transition as a representative of a large class of quasi-one-dimensional materials. 

The theory of spin-Peierls (SP) transition was developed in the 1970's in connection with organic quasi-one- 
dimensional materials. The seminal work of Cross and Fisher^ (CF) is nowadays taken as the canonical theory 
of the spin-Peierls transitions. It is based on a RPA approximation for the phonon field, which corresponds to ne- 
glecting the dynamics of the phonons. Otherwise the low energy spin correlations are exactly taken into account by 
the bosonization method. The CF theory is expected to work when the phonon frequency is much smaller than the 
spin gap. Moreover, CF have shown that corrections to RPA are small if the bare phonon frequency in the direction 
of the magnetic chain is much smaller than the phonon frequency perpendicular to the chains, that is to say if there 
is a preexisting soft mode before coupling phonons with the magnetism. 

In this context the discovery of the first inorganic SP compound CuGeOa opened the possibility to experimentally 
study the dimerized phase and the phase transition with unprecedented precision^. The recent characterization of 
a new family of spin-Peierls materials with an incommensurate intermediate phase has renewed the interest on a 
precise description of the phase transition^. The picture that emerges from neutron scattering measurement of the 
phonon spectra in CuGeOa is the foUowing^'^ : at room temperature, there are two phonon modes relevant for the 
SP transition with wi = 3.12THz and uj2 = 6.53THz at the (0.5 0.5) point of the Brillouin zone. The phonon 
dispersion transverse to the propagation vectors of these relevant phonons is rather flat^. Therefore, no preexisting 
soft mode is observed and the RPA is not a priori justified. Moreover, on lowering the temperature, uji and ^2 do 
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not soften at the SP transition, as expected in the CF scenario. This is a strange behavior for a displacive structural 
phase transition. 

On the other hand, Gros and Werner^ (GW) have identified a parameter regime where the CF approach does not 
predict a softening of the relevant phonons but instead the appearance of a central peak in the dynamical structure 
factor at the transition point. The temperature dependence of the Peierls active phonon modes predicted by this 
analysis are in agreement with neutron scattering determinations in CuGeOs. However the applicability of the 
adiabatic approximation or RPA in this material is not a priori justified. 

In fact, what is missing for a complete understanding of the spin-Peierls transition in CuGeOa, and possibly in 
other inorganic quasi-one-dimensional materials, is a theory which does not rely on the static approximation for the 
phonons. A first step in this direction was presented in [12], where the dynamics of the transverse phonons in the 
antiadiadatic limit was taken into account and the relevance the three dimensional character of the phonons was first 
brought to light. We give some more elements of such a theory in the present work. 

The plan of the paper is as follows. In Section II, we start from a one-dimensional magnetic model coupled to 
two dimensional phonons and integrate out the phonon coordinates. Then we treat the effective in-chain interaction 
exactly by bosonization, and the inter-chain interaction by a mean field approximation. In Section III exact results for 
the ID sine-Gordon modeF are used to solve the mean field equations. By comparing with the adiabatic approximation 
we identify the relevant parameter to analyze the deviations from the adiabatic limit, namely the ratio between the 
frequency of phonon relevant for the SP transition and the frequency of the phonons transverse to the magnetic 
chain {uj±). The adiabatic limit is obtained only for high enough value of Reducing this parameter both the 

gap and the spin-Peierls temperature are reduced. For lu± 0, i.e. the pure one dimensional problem, the gap 
and the SP temperature vanish. We discuss our results in connection with a recent study of the AF chain coupled 
to one-dimensional phonons'^. In Section IV we study the dynamics of the SP relevant phonon mode in the high 
temperature phase, using a RPA on the effective interchain coupling generated by the spin-phonon interaction. We 
then address the question of absence of softening close to the SP transition. Finally, in Section V we present our 
conclusions and prospects for future related work. 



II. THE EFFECTIVE INTERACTIONS INDUCED BY PHONONS 



Let us start by considering a two dimensional system of spin = 1/2 antiferromagnetic Heisenberg chains, with 
exchange constant Jy along one of the axes of a non-deformed square elastic lattice. We label by j chains, and by i 
the sites on each chain. The spin-phonon Hamiltonian is given by^: 

where u\ are the relevant (scalar) coordinates for ion displacements with respect to equilibrium positions, are 
their conjugate momenta, K\\ and K±^ are the elastic couplings along in-chain and inter-chain ions, respectively, and 
a measures the deformation effect on the magnetic exchange constants. 

In the absence of phonons, we can write the bosonized low energy (long distance along the chains) Hamiltonian for 
the Heisenberg spin chains"'^° as 

H^n = • Ci - ^/d^ {kl {d..'^{x))\ {d^^{x)f\ , (2) 

id 

where x = ia (a being the lattice spacing) , (j)^ are the fields dual to the scalar fields 0-' , defined in terms of their 
canonical momentum as dx(t>^ = , and : : stands for normal order with respect to the isolated chains ground 
state. The Fermi velocity of left and right excitations and the usual Luttinger parameter Kl depend on the XX Z 
anisoptropy parameter, taking the values vp = ^J\\'^ ^nd ~ 1 for the isotropic Heisenberg chain with only nearest 
neighbor interaction. A marginal term A : cos{2\/2tt(/){x)) : has been discarded in eq.(2). 

The term coupling spins with phonons can be treated perturbatively by using the continuum expression for the 
spin energy density 

Si ■ ~ p'ix) = -^d.qy'ix) + i-iyp : cos(V2^0^(a;)) : + • • • (3) 

V ZTT 

where /3 is a non-universal constant and the ellipses indicate higher harmonics^^. Notice that : cos(\/27r(/)^) : is a 
relevant operator of scaling dimension D = 1/2. 
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We will also assume that each chain dimcrizcs in such a way that the distortions can be approximately described 
as « (— l)'u-'(x). Notice however that neighboring chains are not correlated a priori. Then 



(4) 



We will keep only the leading order in the gradient expansion. The interaction Hamiltonian, in the continuum limit 
and up to leading order in the lattice spacing, then reads 



Hint = a 



J\\ E«i ~OSl-Sl+, = -gY,jdxu^{x) : cos{V2^<p^ {x)) : , 



(5) 



2a J, ,13 



where g 

Regarding the phonon Hamiltonian, we find it convenient to first construct the corresponding Lagrangian. Up to 
leading order in the lattice spacing the phonon Lagrangian can be written in the continuum limit as 



^f'^ = ^E/'^^ - ^ + \^l){u'{x,t)f + \u:lu^{x,ty+\x,t) 



(6) 



where uof, = 4-^ and lj'^. = 4-^. 

Now we can write the action for the complete system and perform the usual Wick rotation r = it, arriving to the 
Euclidean action describing the low energy physics at finite temperature T: 



S f 
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/dx dr ({drU^X,T)f + ^ + lcul){u^X,T))' - ^UjIu^ + \x , t)u^ {x , t) 

— ffE] / dxdTu^{x,T) : cos(-\/27r0^ (x, r)) :, 
where j runs from 1 to A^, a; from to L, r from to /3 = 1/T, and the spin Euclidean action reads 



Sspin [0"' ] 



dx dT —{dr^'{x,T)y+VFid,^ix,T)y 



(7) 



(8) 



After Fourier transforming j ^ q (the wave vector perpendicular to the magnetic chains) , r — > a;„, with 



the Matsubara frequency, the phonon Green's function reads 

a 1 •^-^ dq e 



2m 2tt uj{q) 



^-^{q)\r~r I ^ 2N{uj{q)) cosh(a;(g)(T' - t)) 



(9) 



where N{uj{q)) is the Bose distribution and the phonon dispersion relation is given by uj'^{q) =0-11+ uj'j_sm{q/2)'^ . 

The action is quadratic in the phonon field and can be easily integrated, rendering the effective action for coupled 
spin chains 

^-// = E Sspinm - y E / / (■ cos(V2^(/)-' (x, r)) : G" (r, r') : cos(\/2^0^" [x, r')) :) . (10) 

3 3,3' 

We find it convenient to separate the in-chain from the inter-chain interactions writing 

•^e// = Sspin[(t>'\ + Sin + Sinter (11) 



where 



Sh 



(12) 
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and 



S^er J J {;■ cos(V2^</)^ (x,t)) : G"'(t,t') : cos{V2^ cjy^ ' {x,t')) :) . 



(13) 



We will now treat the inter-chain interaction with a mean field decoupling. Following Delfino et. al.^^, we write 

: cos{V2^(ty') ■■= eo + N{cos{V2^4)^)), (14) 

where the normal order N is now taken with respect to the vacuum of the coupled chains system. Here eo represents 
the subtraction of the staggered energy density (£0(2; = ia) ^ ( — 1)'+^(< Sf ■ Sl^-^ > — < Sl^-^ ■ Sl^2 >)) ^he 
coupled system, describing the magnetic dimerization along chains. According to the dimerization picture, we take 
eo to be position independent. Then 
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^ y y dT dr' \^el G"' iT,T') + eoN{cos{V2^(j)^ {x,t))) G"' {t,t') 
+ eoiV(cos(y2^0^' (x, r'))) G" (r, t') + • • • } 



(15) 



where the ellipsis corresponds to quadratic fluctuations in the inter-chain coupling, which we will neglect at the MF 
level (i.e. we assume that fluctuations are small with respect to the expectation value eo). Due to isotropy and 
translation symmetry G^^ (ttt') = G'^^^ l(|r — t'\), then the last two terms are equal. Moreover, we can change the 
summation indices from j, j' to j, J = j' — j and r, r' to t, O = r' — r and perform the sum over J and the integral 
over O. We obtain 



r dT'Y,G'{\r-r'\) ^ rrfe[^G^(e)-G°(e)) 



£-{-l[(l-e-'^-«) + 27V(c.||)sinh(/3c.||)] 



I r dq 1 



2 27r Lo{q)^ 
notably leading to the temperature independent result 

r/3 



(1 - e^'^'^(«)) + 2N{uj{q)) sinh(/3cj(g)) }, 



(16) 



/ dr'Y^G'^{\r-r'\) = 
Jo ^0 "i^ll' 



1 - 



(17) 



Then, Sinter reads 



Sinter = "4 y X! J F {ui \\ , UJ ±) 

-eog^Yl / / '^'^ A^(cos(y2^(/)^)) F{uji\,uj±). (18) 
i 

After writing the space-time volume integral, and writing back in the last term N {cos{^/2TT(f>^ (t)) —: cos(-\/27r0^ (r)) : 
—eo, the effective action in eq.(ll) takes the mean field form 

^^f = J^^^P'"^'^^^ - I dx dT dr' : cos(V2^0^(a;,r)) : G°(t,t') : cos(V2^0^(a;, r')) : 

i j 

-tQg^F{u\\,iOx_)Y, j dxdT : cos(V2^0J(x, r)) : +^elg^Fiuj\\,Lu^)NLp. (19) 



The second term in the first line is the remaining in-chain retarded interaction. In the antiadabatic limit (m — s- oo) it 
becomes an instantaneous interaction which is a marginally irrelevant perturbation. Following an argument similar as 
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the one presented in [14] (see the discussion starting from eq.(2.30)), it is possible to show that even in the finite mass 
case this term gives rise to an irrelevant operator (see also [8]). Consistently with the bosonized spin Hamiltonian in 
eq.(2), we neglect it in the following. 
Finally we obtain 

-^e^f = Il5,p,„[0^] - eoff2F(t^||,c^i)^ JdxdT : cos(V2^0^(x, r)) : +el ^^F{u:puj^)NL(3. (20) 

3 J 

The second term is a strongly relevant perturbation that opens a gap for any value of the phonon frequencies t^n and 
uj±, as long as eg does not vanish (this will be evaluated by a self-consistent procedure in the next section). Note 
that this behavior is different than the one obtained in a pure one-dimensional model®. In that case, going from the 
adiabatic to the antiadiabatic limit, the gap closes at an intermediate value of the phonon frequency. Our previous 
discussion shows that when the 2 or 3D character of the phonons is taken into account such a crossover does not show 
up. In [12] a model similar to (1) has been considered only in the antiadiabatic limit; our results show that their 
fully gapped phase extends over a finite frequency range. In addition, as we have kept the explicit dependence on the 
phonon frequencies, we can analyze the evolution of the gap and the SP transition temperature. We will undertake 
this analysis in the following section. 



III. THE SPIN GAP AND TRANSITION TEMPERATURE AS A FUNCTION OF THE PHONON 

FREQUENCY 



Let us start by discussing how the Mean Field adiabatic treatment of the phonon field can be implemented in our 
calculation. This will give results for the spin gap and the transition temperature generalizing those in [13], which are 
in fact extensions of the classical work of Cross and Fisher^. This calculation will set up the energy and temperature 
scale to be compared with our more general non-adiabatic results. Fixing the displacement field in eq.(7) to be a 
constant (chain, position and time independent) u^{x,t) = u, we have 



Sspinifp] ~ / dxdr : cos{V2tt(1>{x,t)) : -\~ — L Pu^i^v^ 



(21) 



Note that the term with ijo±_ vanishes and therefore the total action is N times the action of each chain. 

We obtain it in a self-consistent way by minimizing the total energy with respect to it. The ground state energy 
is evaluated using the exact results for the massive sine-Gordon theory. The lowest energy excitation is a soliton of 
mass Mad, which can be obtained from eq. (12) of [7] { 2fi = gu/vp and /3 = 1/2 in our notation): 



2VF 7rr(|) 



gu 



MadVp 



^r(|)i^ 



4r(i) 



(22) 



The specific ground state energy for the sine-Gordon theory is given by — |;A/^^w|, tan(|-), so that the total energy 
density reads 



ad ^ ^ —uif.u^ - -Mhvliimi-). 

NL 2a II 4 ad F Vgy 



Inverting (22) to obtain Mad{u) and minimizing (23) with respect to u we get: 
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a 1 .TT^TrFd) 

2 2 o ^^^^7^) 



Avp mwi 3 




H 



g^_a_ 
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where H — 5.4133 is just a constant collecting all of the numerical factors. The spin gap A"'' is then given by 



(23) 



(24) 



(25) 



Equations (24) and (25) give the zero temperature gap in terms of the model parameters g, J\\ and cij|| . It is equivalent 
to the first of the eqs. (8) in^^. 



6 



In the Mean Field adiabatic treatment, the spin-Peierls transition temperature Tgp can be obtained by considering 
the finite temperature free energy. Borrowing the appropriate expansion in powers of u from Eq. (12) of [13] {j^^ 
must be replaced by gu/ ^fa) the lowest order of the sine-Gordon model free energy density reads 

so that the full variational free energy is obtained by adding the contribution of the elastic energy, 



pad 771 TT 9 TT'^ (quY , As , 



The transition temperature corresponds to the one where the free energy minimum shifts from u ^ to u ~ 0. 
Although we have not written quartic terms in u, it is easy to see that the transition is signaled by a change of sign 
in the coefficient of u^. This condition gives 

2 2 

We now undertake the calculation of the gap and the transition temperature in the non adiabatic case. We start 
from the action (20), our main result previously obtained. Comparing cqs. (20) and (21) we note that the soliton mass 
AI can be computed by following similar steps as we did to obtain Af"''. Now, the mean field variational parameter 
is eg instead of u (indeed the factor gFe^ plays the role of u); the total energy density is given by 

' " = - ^M^4tan(^). (29) 

We readily obtain the soliton mass as 

M = H^F{uji\,uji_), (30) 

and the corresponding gap A = Mvp. Comparing this result with eq.(24) we obtain the relation between the gap in 
the non-adiabatic and the adiabatic case: 

-^^1- , ' . (31) 



In Fig. 1 we show the evolution of the gap with the ratio Only for ^ — > cxd, i.e. when the lines perpendicular to 

the chain move rigidly, the gap corresponds to the adiabatic one. Note that this is precisely the condition established 
by CF for the validity of their RPA calculations. In the general case the gap decreases when ^ decreases. This is 
an important result applicable to a general non-adiabatic spin-Peierls system. For example, for CuGeOs we can take 
as the relevant phonon mode (with frequency ix>||) the projections of the (0.5 0.5) phonons on the direction of the 
magnetic chains. The transverse direction is then the ^-direction along the path (0.5 x 0.5). The phonon dispersion 
relation uj'^{q) = + uj'j_ sin(g/2)^ in our model corresponds to phonons in this direction. We estimate uj± from the 
frequencies at the center and at the boundary of the Brillouin zone in this direction. It has been recognized that 
there are two phonons at (0.5 0.5) which are the relevant for the SP distortion (with frequencies cui = 3.12THz 
and UJ2 = 6.53 Ti/z at room temperature)^. From Fig. 12 in [5] we can follow the dispersion of the branches in the 
direction (0.5 x 0.5) connected with these two modes. The lowest energy branch is rather flat and we do not take 
it into account. For the phonon with the frequency UJ2 we have ^ 0.58 and the gap from eq. (31) is reduced 

to 13% of the one estimated by the adiabatic approximation. In practice, the measured zero temperature gap is 
used to estimate the spin-phonon coupling^^. Our previous result implies that the spin-phonon coupling could be 
strongly underestimated by using adiabatic results. Indeed, from our present approach we estimate the dimensionless 
spin-phonon coupling (A = ^^^) for CuGcOa one order of magnitude larger than the previously calculated one. This 
could have important consequences to interpret experimental results in this and other materials. 

A relation similar to eq. (31) holds between the transition temperature in the non-adiabatic and the adiabatic case. 
The free energy corresponding to action (20) should be expanded in powers of cq. We take again this expansion from 
Eq. (12) of [13] (now ^^'^^^2 is replaced by egg'^ F{ll>\^, ujj_) / y/a). The total free energy density is 



.J^^ _J^rp2 , 2 g'-P^(^lh^^) A _ .9^7r^f(^||,^x) \ 
■^^ NL 6vF ° 2 1^ 2r4(|)T J 



(32) 
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FIG. 1: Dependence of the gap and the SP temperature with 



The vanishing of the factor proportional to eg signals the transition temperature. We get 



and the promised relation is: 



9 9 

J^SP ^ F{uj\\,uj±_) 



Tsp 

rpad 



(33) 



(34) 



Notice that the ratio from our computation docs not depend on The value of such a ratio seems to be 
fixed by the adiabatic calculation. The existence of such universal ratio is similar to the BCS mean field theory of 
superconductivity. Moreover, the experimental value observed in CuGeOa is close to the BCS-ratio 1.76. In fact 
a BCS-like theory could be obtained from the original Heisenberg model coupled to phonons by a Jordan- Wigner 
transformation and neglecting interaction between the spinless fermions. Note however that the spinless fermions are 
strongly interacting, and this fact renormalizes the BCS ratio away from the non interacting value^'^. The precise value 
of this ratio is difficult to obtain because it is affected by logarithmic corrections induced by marginally irrelevant 
terms. Indeed, the effect of the marginal term arising from the Heisenberg interaction is expected to be small in 
CuGeOa because the second next- nearest neighbor is closer to the critical value^^, where logarithmic corrections 
vanish. However, as can be seem from action (11), logarithmic corrections are also expected from the spin-phonon 
interaction term and they do depend on the phonon frequencies. Therefore we expect that the exact value of -jA^ will 
ultimately depend on — . 



IV. DYNAMICS OF THE PHONON MODES 



The original work of CF on spin-Pcierls systems predicts the softening of a phonon mode above the transition 
temperature, whose polarization pattern corresponds to the static distortion leading to the low temperature phase 
below the transition point Tsp. The softening is due to the coupling of the phonons with the one-dimensional 
magnetism. This is in fact a quite general scenario for a structural displacive phase transition. Therefore it came as 
a surprise that the SP transition in CuGeOa shows no phonon softening^^'^°. Gros and Werner^ have reanalyzed the 
CF approach showing that, for higher enough bare phonon frequency, there is a second pole of the phonon propagator 
without softening of the active SP phonon, characterized by the presence of a central peak in the dynamical structure 
factor near the transition point. Moreover, the RPA approach underlying the CF calculation is justified if the 
SP phonon has a bare frequency considerably lower than other momentum phonons, i.e. a non-magnetic softening 
observable at high temperature relative to the SP transition temperature. This is not the situation in CuGeOa and 
application of RPA seems questionable. 

In the following we show how by extending our previous calculations to the dynamics of the phonon modes we 
can solve this puzzle. The extension of the static MF approximation on the interchain interaction to the dynamical 
correlation function is accomplished by a RPA analysis along the lines of [12.21]. In contrast to the one-dimensional 



8 



analysis in CF. this RPA is based on the gMasi-one-dimensionahty underlying this kind of materials, and does not 
require a preexisting mode softening. It becomes exact in the limit of an infinite number of neighboring chains. 

The dressed phonon frequencies are obtained from the poles of the retarded phonon Green's function D{k,q,ujy^ 
which can be computed from the Matsubara Green's function D{k, g, a;„), with wave vectors k in the chain direction 
and q the transverse one, by analytic continuation as D{k,q,uj) = D{k, q^icon uj). This is in fact the Fourier 
transform of the Green's function in real space and imaginary time D^~^ {x — x',r — r'), which can be found by 
functional derivation of the partition function when a current term is added: 



Z[J]= J 2?u^I?(/)^exp|-^ J dxdT[C{(fy>{x,T),u^{x,T))+gJ^{x,T)u^{x,T)]^ 
Integrating out the phonon coordinates as before we obtain the following effective action 



(35) 



2 



3, J 

The dressed phonon propagator 

(37) 



{JHx.r)}=0 



g2 SJi'{x',T')SJi{x,T) 
can be obtained from eqs. (35, 36, 37) as 

D^-''{x-x',t-t') = G^'-^{x' -x,t' -T)-g^Y2 J dxdx'drdT'G'^-' {x - x,t - t)g'^' {i' - x' ,? - t') 

hi' 

(t(: cos{V2^(j^{x,T)) :: cos{V2^(j/' {x' ,?)) :)), (38) 



where G^-J {x - x' ,t - t') ^ 6{x - x')G" (t, t') is a convenient notation for the non-interacting phonon propagator 
and T indicates Euclidean time-ordering. After Fourier transforming to [k^q,ojn) coordinates and performing the 
analytic continuation iw„ — > w we obtain the desired retarded phonon Green's function 

D{k,q,Lo) = G{q,Lu)+g^G\q,Lu)x{k,q,uj), (39) 

where G{q, cu) = -~uj'^+ui(q) ^'^'^ bare phonon Green's function and x(^j ^) stems for the Fourier transform of the 



retarded correlator of : cos(v27r(/)^) : . In order to evaluate bulk corrections arising from inter-chain coupling, we will 
calculate this correlator using a RPA approach in the interchain interaction^^. It reads 

x^^^l^-?,-) - 1 ^r^'^^t^'fl, rv (40) 

1 - g^Ginter[q,uj)x"[k,uj;l ) 

where x°(fc)'i';T') is the already known finite temperature one chain : cos{^/2tt(1>) :-correlator and Ginter is the inter- 
chain bare phonon propagator, defined in coordinate space by the subtraction 

GlnLi^ - x', T - r') = G^~i'{x ~x',T- t') - P^'G°{x ~x',t- t'). (41) 

The one chain : cos(-\/27r0) :-correlator at finite temperature has been computed elsewhere^'^'^^, being given by 

where d ~ 0.37 is a constant weakly depending of the momentum cutoff in the bosonization procedure, A = vp\k — 7r| 
is the lower edge of the two spinon continuum, and 



1 r 

Jo 



dxe^''^ {smh{x)) ^ 



1 Ffi 



4 " 2 



2x/2^r[|-i|] 



(43) 
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It will be useful to note that at low frequency one can expand 

Tx°(fc,.;T)«.xo-.Xi^+X.(^)\ (44) 

with xo 0.26 and X2 ~ 2.2. 

The Fourier transform of the interchain propagator G*"*'^'' is obtained (see Section II) as 

G-*-(g,c.) = ,\ . -T{uj), (45) 



with 



2tt 



We are interested in the temperature evolution of the SP active phonon mode, given by fc = tt and q = 0. From 
eq. (39-45) we have 

D(0, .) = ..-(l-..-T(.)x°(.,c.;r)) 

G-i + ^(l-G-iT(c.))x''(^,u;;T) 

and the phonon frequencies are given by the poles over the real axis of this expression, i.e. the roots of the following 
equation 

^ ^2 ^ _ (^2 _ iu^)T{u;))Rex°{7:, co; T) = 0. (48) 

Let us start the analysis of the consequences of the previous results by recalculating the spin-Peierls transition 
temperature. Note that, different than the previous Section, we are now coming from the high temperature phase to 
the transition point. As discussed by GW, Tgp is signaled by a macroscopic occupation of the Peierls active phonon 
mode, i.e. the transition temperature takes place when cq. (48) has a solution lj = Q. We then obtain 



^^^^ 



1 - 



(49) 



It is important to stress that eq. (49) gives the same dependency of Tsp on the microscopic parameters than eq. (33) 
in the previous Section; differences in the numerical prefactor arise from different criteria to fix the momentum cutoff. 
Note however that relation (34) is fulfilled. This results prove that our RPA procedure contains the previous MF 
approach when the static properties are considered. 

In the rest of this Section we face the question of phonon softening, discussing the calculation of the renormalized 
phonon frequencies. We use (49) to eliminate the cut-off dependent coupling from cq. (48). After straightforward 
algebra we obtain 



XO 1- , - -1 / 1- . ^ " Ur,pfiexV,c.). (50) 




This equation is the generalization of the GW results to finite values of and agrees with their corresponding eq. 

(6) when ^ ^ 1. There are two different temperature dependence regimes of the renormalized phonon frequency 

^(TrjO), as depicted in Fig. 2. For low enough ti>|| the phonon progressively softens from higher T until Tsp, where it 
drops to zero frequency. In this regime the phonon softening explains the SP transition. For higher W|| , the phonon 
frequency remains finite up to the SP temperature and no soft phonon behavior is obtained; in fact, there can be 
partial softening or even hardening of the renormalized frequency close to Tsp. However, an additional solution lo = 
in eq. (50) signals the appearance of a central peak in the dynamical structure factor, leading to the SP transition^. 

The limiting value of uj\\/Tsp for separation of regimes depends on lo^. To find this transition analytically it is 
enough to compare concavities of both sides of eq. (50) at T = Tsp, as seen in Fig. 2. Expanding both sides up 
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FIG. 2: Evolution of the SP relevant phonon mode. Dashed lines, from blue to red, represent the rhs of eq. (50) when 
temperature is lowered until Tsp. Solid gray lines represent the Ihs for lj^^/Tsp — 2,3,4 and uij_/Tsp = 1. The intersections 
illustrate the evolution of the renormalized phonon frequency. For ui^^/Tsp — 2 the phonon softens completely, for ui^^/Tsp = 3 
it softens partially, and for ijj\\/Tsp = 4 it hardens close to Tsp- 



to second order in the soft phonon regime is reached when, at T = Tsp, the coefficient of ijP in the r.h.s. is 

smaller than the one in the l.h.s. Under this condition the bare phonon frequency (that obtained at hight T) evolves 
continuously to zero when the temperature is lowered until the transition temperature. The following inequality 
should then be fulfilled to have the soft phonon regime: 



-11 iUi^^TT^U^^ ^^^^ 



Tsp V X2 V 2 1 + {uj±/u:\\Y 



In the other case, the bare phonon frequency evolves to some finite value at the SP transition. In Fig. 3 we show 
the phase diagram of the spin-phonon system on the uj±^ — lo\\ plane. Note that in a wide range of high enough 

the separatrix is given by ^ 2.2 as in GW. Our result extends this condition for finite values of where 
realistic spin-Peierls materials as tJuGeOa live. This also explains why the RPA approach for the spin-phonon coupling 
underlying the GW calculation compares well with experiments on CuGeOa. 

In Fig. 4 we show the temperature dependence of the renormalized phonon frequency for ~ 1 where the 
phase diagram of Fig. 3 is not completely flat. The behavior shifts from a soft to a non soft phonon regime at some 
intermediate value of ;f^- In this last regime the mechanism of the SP transition has been associated with the 
emergence of a central peak in the spectral signal measured in neutron scatcring experiments^^. For high enough 

the phonon hardens from a minimum value at some intermediate temperature. kX, T = Tgp the frequency of 
the phonon is the bare one as in the high temperature limit. We will use this fact to fix uj^ for CuGeOa. Note that 
the T-dependence of the phonon is similar to the one obtained by GW but the frequency and temperature scale are 
different. 

From the previous results we can follow the T-dependence of the Peierls active phonons in CuGeOa. In Fig. 4 (6), 
we show our prediction for the higher energy mode T^, which is the most anomalous and strongly coupled to the 
magnetism"*. We take as a bare phonon frequency lo\\ ~ 6.8THz i.e. its low temperature value. For uj± we use the 
relationship obtained in the previous Section. Note that a value of Tsp = 14.1ii'(0.294TiJz) has been used to give 
the frequencies in Thz and the temperatures in K. Our results compare well with Fig. 3 of [4]. In fact we predict a 
stronger hardening than the one experimentally observed (7% against 4.5% seen in the experiment). This could be 
due to the simplified model we use, which does not take into account the fact that two phonon modes are necessary 
to describe the SP transition in CuGcOa"'. 

Besides, we can make a prediction for the phonon frequency behavior at higher temperature than the one already 
measured. It should increase as shown in Fig. 4, approaching the bare frequency at very high T. We note that at 
room temperature the phonon frequency does not correspond to the bare one as expected. It is almost the minimum 
of the predicted phonon frequency. The bare frequency is only obtained at T ^ Tsp and at very high temperature 
(probably not experimentally accessible). 

To summarize the results of the present Section, we have generalized the MF approach of the previous Section to 
the dynamical correlation functions in the high temperature phase. Our approach relies on a RPA on the effective 
interchain coupling generated by the spin-phonon interaction. In contrast to the previous theoretical calculations. 
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No Soft Phonon 




Soft Phonon 



0.5 - 



FIG. 3; The line separating the zone in the t^^, plane, where the spin-Peierls transition takes place by softening of the 
phonon from the one where the Peierls active phonon does not soften completely. 



03, =3.94 THz m,,=6.8 THz 




FIG. 4: (a)The temperature dependence of the renormalized phonon frequency for different values of the bare with 

(b) The hardening of the phonon for the parameters describing the higher energy active SP mode in CuGeOs. The bare values 

used for clI|| and cjx are shown on top of the Figure. 



the present one does not presuppose the existence of a soft phonon mode in the non magnetic regime. Therefore, it 
can been apphed to CuGeOa. Moreover, the quaUtative behavior of the Peierls active phonon is similar to the one 
previously obtained. There is a regime corresponding to low frequency phonon where the structural phase transition 
takes place by progressive softening of a phonon from the high temperature phase. For higher bare phonon frequency 
no soft phonon is observed and it can even harden from the room temperature value to the transition temperature 
one. This is precisely the behavior seen in neutron scattering measurements of the phonon spectra in CuGeOa. 



V. CONCLUSIONS 



We have developed a route for the spin-Peierls transition which goes beyond the usual adiabatic treatment of 
the phonon field. We emphasized the essential character that plays the dispersion of the phonon in the transversal 
direction to the magnetic chains. Moreover, the results of the present paper show that, when the adiabatic hypothesis 
for the phonon coordinates is relaxed, the one-chain model does not represent a good starting point to describe a 
real system. Furthermore the effective interchain interaction generated by the phonons is an essential ingredient for 
the opening of the gap and the existence of a finite temperature phase transition. The precise determination of the 
SP transition temperature from the microscopic parameters strongly depends on the width of the transversal phonon 
dispersion. 

Our approach for the dynamics of the phonon mode in the high temperature phase justifies and generalizes previous 
calculations of the phonon spectra based on the Cross and Fisher early work. Applied to CuGeOa, we show that 
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the phonon does not soften until the phase transition. Moreover the overall temperature dependence of the phonon 
frequency is consistent with experimental determinations. 

As the elastic interaction between the chains play a central rol in non-adiabatic spin-Peierls systems, we anticipate 
that a change in the relative positions of the magnetic ions in different chains will change the phase diagram and 
the mechanism of the phase transition. In this sense it is highly interesting to study the recently discovered quasi- 
one-dimensional magneto-elastic system TiOCl with an intermediate incommensurate phase where the Ti ions are 
not aligned in the direction perpendicular to the magnetic chain'^. We will study in a forthcoming work this system 
within the formalism developed in the present paper. 
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We develop a theory of the spin-Peierls transition taking into account the three dimensional 
character of the phonon field. Our approach does not rely on the adiabatic or mean field treatment 
for the phonons. It is instead based in the exact integration of the phonon field, the exact long 
wavelength solution of the one chain spin problem, and then a mean field approximation for the 
interchain interaction. We show that the spin gap and the critical temperature are strongly reduced 
due to the finite frequency effects of the phonon coupling transverse to the magnetic chains. We 
claim that our results should be applicable to the inorganic spin-Peierls compound CuGeOs. We 
show that the long standing discussion on absence of a soft mode in this compound can be naturally 
resolved within our theory. 

PACS numbers: 75.10.Dg, 75.10.Jm, 63.70.+h 

I. INTRODUCTION 

Quasi-one-dimensional electronic and magnetic materials are fascinating objects which have been intensively studied 
during the last decades. They show a variety of different phases as the magnetically dimerized, unconventional metal 
and triplet superconductor phases. It is possible to go from one phase to another by changing the experimentally 
accessible macroscopic parameters: pressure, temperature and the applied magnetic field. The different phases arise 
from the competition between the Coulomb repulsion between electrons and their coupling with the lattice degrees of 
freedom. Whereas the essential physics relics on the one-dimensionality of the electronic dynamics, the stabilization of 
the different phases depends on the three dimensional character of the system and thus on the interchain interaction. 
Transitions between the different phases are, in a way, a dimensional crossover from ID to the three dimensionality 
of any real material. 

Spin-Peierls systems are probably the simplest example of quasi-one-dimensional materials. They show a high 
temperature phase which can be well modelled by a one dimensional antiferromagnetic Heisenberg chain and a low 
temperature phase with lattice dimerization and spin singlet ground state. The transition between them arises from 
the competition between the gain in antiferromagnetic energy and the loss in elastic energy. It is indeed possible 
due to the three dimensional character of the phonon field. It is therefore very attractive to carefully understand the 
precise character of this transition as a representative of a large class of quasi-one-dimensional materials. 

The theory of spin-Peierls (SP) transition was developed in the 1970's in connection with organic quasi-one- 
dimensional materials. The seminal work of Cross and Fisher^ (CF) is nowadays taken as the canonical theory 
of the spin-Peierls transitions. It is based on a RPA approximation for the phonon field, which corresponds to ne- 
glecting the dynamics of the phonons. Otherwise the low energy spin correlations are exactly taken into account by 
the bosonization method. The CF theory is expected to work when the phonon frequency is much smaller than the 
spin gap. Moreover, CF have shown that corrections to RPA are small if the bare phonon frequency in the direction 
of the magnetic chain is much smaller than the phonon frequency perpendicular to the chains, that is to say if there 
is a preexisting soft mode before coupling phonons with the magnetism. 

In this context the discovery of the first inorganic SP compound CuGeOa opened the possibility to experimentally 
study the dimerized phase and the phase transition with unprecedented precision^. The recent characterization of 
a new family of spin-Peierls materials with an incommensurate intermediate phase has renewed the interest on a 
precise description of the phase transition^. The picture that emerges from neutron scattering measurement of the 
phonon spectra in CuGeOa is the foUowing^'^ : at room temperature, there are two phonon modes relevant for the 
SP transition with wi = 3.12THz and uj2 = 6.53THz at the (0.5 0.5) point of the Brillouin zone. The phonon 
dispersion transverse to the propagation vectors of these relevant phonons is rather flat^. Therefore, no preexisting 
soft mode is observed and the RPA is not a priori justified. Moreover, on lowering the temperature, uji and ^2 do 



2 



not soften at the SP transition, as expected in the CF scenario. This is a strange behavior for a displacive structural 
phase transition. 

On the other hand, Gros and Werner^ (GW) have identified a parameter regime where the CF approach does not 
predict a softening of the relevant phonons but instead the appearance of a central peak in the dynamical structure 
factor at the transition point. The temperature dependence of the Peierls active phonon modes predicted by this 
analysis are in agreement with neutron scattering determinations in CuGeOs. However the applicability of the 
adiabatic approximation or RPA in this material is not a priori justified. 

In fact, what is missing for a complete understanding of the spin-Peierls transition in CuGeOa, and possibly in 
other inorganic quasi-one-dimensional materials, is a theory which does not rely on the static approximation for the 
phonons. A first step in this direction was presented in [12], where the dynamics of the transverse phonons in the 
antiadiadatic limit was taken into account and the relevance the three dimensional character of the phonons was first 
brought to light. We give some more elements of such a theory in the present work. 

The plan of the paper is as follows. In Section II, we start from a one-dimensional magnetic model coupled to 
two dimensional phonons and integrate out the phonon coordinates. Then we treat the effective in-chain interaction 
exactly by bosonization, and the inter-chain interaction by a mean field approximation. In Section III exact results for 
the ID sine-Gordon modeF are used to solve the mean field equations. By comparing with the adiabatic approximation 
we identify the relevant parameter to analyze the deviations from the adiabatic limit, namely the ratio between the 
frequency of phonon relevant for the SP transition and the frequency of the phonons transverse to the magnetic 
chain {uj±). The adiabatic limit is obtained only for high enough value of Reducing this parameter both the 

gap and the spin-Peierls temperature are reduced. For lu± 0, i.e. the pure one dimensional problem, the gap 
and the SP temperature vanish. We discuss our results in connection with a recent study of the AF chain coupled 
to one-dimensional phonons'^. In Section IV we study the dynamics of the SP relevant phonon mode in the high 
temperature phase, using a RPA on the effective interchain coupling generated by the spin-phonon interaction. We 
then address the question of absence of softening close to the SP transition. Finally, in Section V we present our 
conclusions and prospects for future related work. 



II. THE EFFECTIVE INTERACTIONS INDUCED BY PHONONS 



Let us start by considering a two dimensional system of spin = 1/2 antiferromagnetic Heisenberg chains, with 
exchange constant Jy along one of the axes of a non-deformed square elastic lattice. We label by j chains, and by i 
the sites on each chain. The spin-phonon Hamiltonian is given by^: 

where u\ are the relevant (scalar) coordinates for ion displacements with respect to equilibrium positions, are 
their conjugate momenta, K\\ and K±^ are the elastic couplings along in-chain and inter-chain ions, respectively, and 
a measures the deformation effect on the magnetic exchange constants. 

In the absence of phonons, we can write the bosonized low energy (long distance along the chains) Hamiltonian for 
the Heisenberg spin chains"'^° as 

H^n = • Ci - ^/d^ {kl {d..'^{x))\ {d^^{x)f\ , (2) 

id 

where x = ia (a being the lattice spacing) , (j)^ are the fields dual to the scalar fields 0-' , defined in terms of their 
canonical momentum as dx(t>^ = , and : : stands for normal order with respect to the isolated chains ground 
state. The Fermi velocity of left and right excitations and the usual Luttinger parameter Kl depend on the XX Z 
anisoptropy parameter, taking the values vp = ^J\\'^ ^nd ~ 1 for the isotropic Heisenberg chain with only nearest 
neighbor interaction. A marginal term A : cos{2\/2tt(/){x)) : has been discarded in eq.(2). 

The term coupling spins with phonons can be treated perturbatively by using the continuum expression for the 
spin energy density 

Si ■ ~ p'ix) = -^d.qy'ix) + i-iyp : cos(V2^0^(a;)) : + • • • (3) 

V ZTT 

where /3 is a non-universal constant and the ellipses indicate higher harmonics^^. Notice that : cos(\/27r(/)^) : is a 
relevant operator of scaling dimension D = 1/2. 
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We will also assume that each chain dimcrizcs in such a way that the distortions can be approximately described 
as « (— l)'u-'(x). Notice however that neighboring chains are not correlated a priori. Then 



(4) 



We will keep only the leading order in the gradient expansion. The interaction Hamiltonian, in the continuum limit 
and up to leading order in the lattice spacing, then reads 



Hint = a 



J\\ E«i ~OSl-Sl+, = -gY,jdxu^{x) : cos{V2^<p^ {x)) : , 



(5) 



2a J, ,13 



where g 

Regarding the phonon Hamiltonian, we find it convenient to first construct the corresponding Lagrangian. Up to 
leading order in the lattice spacing the phonon Lagrangian can be written in the continuum limit as 



^f'^ = ^E/'^^ - ^ + \^l){u'{x,t)f + \u:lu^{x,ty+\x,t) 



(6) 



where uof, = 4-^ and lj'^. = 4-^. 

Now we can write the action for the complete system and perform the usual Wick rotation r = it, arriving to the 
Euclidean action describing the low energy physics at finite temperature T: 



S f 



3 

m 
2 



/dx dr ({drU^X,T)f + ^ + lcul){u^X,T))' - ^UjIu^ + \x , t)u^ {x , t) 

— ffE] / dxdTu^{x,T) : cos(-\/27r0^ (x, r)) :, 
where j runs from 1 to A^, a; from to L, r from to /3 = 1/T, and the spin Euclidean action reads 



Sspin [0"' ] 



dx dT —{dr^'{x,T)y+VFid,^ix,T)y 



(7) 



(8) 



After Fourier transforming j ^ q (the wave vector perpendicular to the magnetic chains) , r — > a;„, with 



the Matsubara frequency, the phonon Green's function reads 

a 1 •^-^ dq e 



2m 2tt uj{q) 



^-^{q)\r~r I ^ 2N{uj{q)) cosh(a;(g)(T' - t)) 



(9) 



where N{uj{q)) is the Bose distribution and the phonon dispersion relation is given by uj'^{q) =0-11+ uj'j_sm{q/2)'^ . 

The action is quadratic in the phonon field and can be easily integrated, rendering the effective action for coupled 
spin chains 

^-// = E Sspinm - y E / / (■ cos(V2^(/)-' (x, r)) : G" (r, r') : cos(\/2^0^" [x, r')) :) . (10) 

3 3,3' 

We find it convenient to separate the in-chain from the inter-chain interactions writing 

•^e// = Sspin[(t>'\ + Sin + Sinter (11) 



where 



Sh 



(12) 
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and 



S^er J J {;■ cos(V2^</)^ (x,t)) : G"'(t,t') : cos{V2^ cjy^ ' {x,t')) :) . 



(13) 



We will now treat the inter-chain interaction with a mean field decoupling. Following Delfino et. al.^^, we write 

: cos{V2^(ty') ■■= eo + N{cos{V2^4)^)), (14) 

where the normal order N is now taken with respect to the vacuum of the coupled chains system. Here eo represents 
the subtraction of the staggered energy density (£0(2; = ia) ^ ( — 1)'+^(< Sf ■ Sl^-^ > — < Sl^-^ ■ Sl^2 >)) ^he 
coupled system, describing the magnetic dimerization along chains. According to the dimerization picture, we take 
eo to be position independent. Then 



' 2 



^ y y dT dr' \^el G"' iT,T') + eoN{cos{V2^(j)^ {x,t))) G"' {t,t') 
+ eoiV(cos(y2^0^' (x, r'))) G" (r, t') + • • • } 



(15) 



where the ellipsis corresponds to quadratic fluctuations in the inter-chain coupling, which we will neglect at the MF 
level (i.e. we assume that fluctuations are small with respect to the expectation value eo). Due to isotropy and 
translation symmetry G^^ (ttt') = G'^^^ l(|r — t'\), then the last two terms are equal. Moreover, we can change the 
summation indices from j, j' to j, J = j' — j and r, r' to t, O = r' — r and perform the sum over J and the integral 
over O. We obtain 



r dT'Y,G'{\r-r'\) ^ rrfe[^G^(e)-G°(e)) 



£-{-l[(l-e-'^-«) + 27V(c.||)sinh(/3c.||)] 



I r dq 1 



2 27r Lo{q)^ 
notably leading to the temperature independent result 

r/3 



(1 - e^'^'^(«)) + 2N{uj{q)) sinh(/3cj(g)) }, 



(16) 



/ dr'Y^G'^{\r-r'\) = 
Jo ^0 "i^ll' 



1 - 



(17) 



Then, Sinter reads 



Sinter = "4 y X! J F {ui \\ , UJ ±) 

-eog^Yl / / '^'^ A^(cos(y2^(/)^)) F{uji\,uj±). (18) 
i 

After writing the space-time volume integral, and writing back in the last term N {cos{^/2TT(f>^ (t)) —: cos(-\/27r0^ (r)) : 
—eo, the effective action in eq.(ll) takes the mean field form 

^^f = J^^^P'"^'^^^ - I dx dT dr' : cos(V2^0^(a;,r)) : G°(t,t') : cos(V2^0^(a;, r')) : 

i j 

-tQg^F{u\\,iOx_)Y, j dxdT : cos(V2^0J(x, r)) : +^elg^Fiuj\\,Lu^)NLp. (19) 



The second term in the first line is the remaining in-chain retarded interaction. In the antiadabatic limit (m — s- oo) it 
becomes an instantaneous interaction which is a marginally irrelevant perturbation. Following an argument similar as 
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the one presented in [14] (see the discussion starting from eq.(2.30)), it is possible to show that even in the finite mass 
case this term gives rise to an irrelevant operator (see also [8]). Consistently with the bosonized spin Hamiltonian in 
eq.(2), we neglect it in the following. 
Finally we obtain 

-^e^f = Il5,p,„[0^] - eoff2F(t^||,c^i)^ JdxdT : cos(V2^0^(x, r)) : +el ^^F{u:puj^)NL(3. (20) 

3 J 

The second term is a strongly relevant perturbation that opens a gap for any value of the phonon frequencies t^n and 
uj±, as long as eg does not vanish (this will be evaluated by a self-consistent procedure in the next section). Note 
that this behavior is different than the one obtained in a pure one-dimensional model®. In that case, going from the 
adiabatic to the antiadiabatic limit, the gap closes at an intermediate value of the phonon frequency. Our previous 
discussion shows that when the 2 or 3D character of the phonons is taken into account such a crossover does not show 
up. In [12] a model similar to (1) has been considered only in the antiadiabatic limit; our results show that their 
fully gapped phase extends over a finite frequency range. In addition, as we have kept the explicit dependence on the 
phonon frequencies, we can analyze the evolution of the gap and the SP transition temperature. We will undertake 
this analysis in the following section. 



III. THE SPIN GAP AND TRANSITION TEMPERATURE AS A FUNCTION OF THE PHONON 

FREQUENCY 



Let us start by discussing how the Mean Field adiabatic treatment of the phonon field can be implemented in our 
calculation. This will give results for the spin gap and the transition temperature generalizing those in [13], which are 
in fact extensions of the classical work of Cross and Fisher^. This calculation will set up the energy and temperature 
scale to be compared with our more general non-adiabatic results. Fixing the displacement field in eq.(7) to be a 
constant (chain, position and time independent) u^{x,t) = u, we have 



Sspinifp] ~ / dxdr : cos{V2tt(1>{x,t)) : -\~ — L Pu^i^v^ 



(21) 



Note that the term with ijo±_ vanishes and therefore the total action is N times the action of each chain. 

We obtain it in a self-consistent way by minimizing the total energy with respect to it. The ground state energy 
is evaluated using the exact results for the massive sine-Gordon theory. The lowest energy excitation is a soliton of 
mass Mad, which can be obtained from eq. (12) of [7] { 2fi = gu/vp and /3 = 1/2 in our notation): 



2VF 7rr(|) 



gu 



MadVp 



^r(|)i^ 



4r(i) 



(22) 



The specific ground state energy for the sine-Gordon theory is given by — |;A/^^w|, tan(|-), so that the total energy 
density reads 



ad ^ ^ —uif.u^ - -Mhvliimi-). 

NL 2a II 4 ad F Vgy 



Inverting (22) to obtain Mad{u) and minimizing (23) with respect to u we get: 
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a 1 .TT^TrFd) 

2 2 o ^^^^7^) 



Avp mwi 3 




H 



g^_a_ 

9 9 

V% mLuf, 



where H — 5.4133 is just a constant collecting all of the numerical factors. The spin gap A"'' is then given by 



(23) 



(24) 



(25) 



Equations (24) and (25) give the zero temperature gap in terms of the model parameters g, J\\ and cij|| . It is equivalent 
to the first of the eqs. (8) in^^. 
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In the Mean Field adiabatic treatment, the spin-Peierls transition temperature Tgp can be obtained by considering 
the finite temperature free energy. Borrowing the appropriate expansion in powers of u from Eq. (12) of [13] {j^^ 
must be replaced by gu/ ^fa) the lowest order of the sine-Gordon model free energy density reads 

so that the full variational free energy is obtained by adding the contribution of the elastic energy, 



pad 771 TT 9 TT'^ (quY , As , 



The transition temperature corresponds to the one where the free energy minimum shifts from u ^ to u ~ 0. 
Although we have not written quartic terms in u, it is easy to see that the transition is signaled by a change of sign 
in the coefficient of u^. This condition gives 

2 2 

We now undertake the calculation of the gap and the transition temperature in the non adiabatic case. We start 
from the action (20), our main result previously obtained. Comparing cqs. (20) and (21) we note that the soliton mass 
AI can be computed by following similar steps as we did to obtain Af"''. Now, the mean field variational parameter 
is eg instead of u (indeed the factor gFe^ plays the role of u); the total energy density is given by 

' " = - ^M^4tan(^). (29) 

We readily obtain the soliton mass as 

M = H^F{uji\,uji_), (30) 

and the corresponding gap A = Mvp. Comparing this result with eq.(24) we obtain the relation between the gap in 
the non-adiabatic and the adiabatic case: 

-^^1- , ' . (31) 



In Fig. 1 we show the evolution of the gap with the ratio Only for ^ — > cxd, i.e. when the lines perpendicular to 

the chain move rigidly, the gap corresponds to the adiabatic one. Note that this is precisely the condition established 
by CF for the validity of their RPA calculations. In the general case the gap decreases when ^ decreases. This is 
an important result applicable to a general non-adiabatic spin-Peierls system. For example, for CuGeOs we can take 
as the relevant phonon mode (with frequency ix>||) the projections of the (0.5 0.5) phonons on the direction of the 
magnetic chains. The transverse direction is then the ^-direction along the path (0.5 x 0.5). The phonon dispersion 
relation uj'^{q) = + uj'j_ sin(g/2)^ in our model corresponds to phonons in this direction. We estimate uj± from the 
frequencies at the center and at the boundary of the Brillouin zone in this direction. It has been recognized that 
there are two phonons at (0.5 0.5) which are the relevant for the SP distortion (with frequencies cui = 3.12THz 
and UJ2 = 6.53 Ti/z at room temperature)^. From Fig. 12 in [5] we can follow the dispersion of the branches in the 
direction (0.5 x 0.5) connected with these two modes. The lowest energy branch is rather flat and we do not take 
it into account. For the phonon with the frequency UJ2 we have ^ 0.58 and the gap from eq. (31) is reduced 

to 13% of the one estimated by the adiabatic approximation. In practice, the measured zero temperature gap is 
used to estimate the spin-phonon coupling^^. Our previous result implies that the spin-phonon coupling could be 
strongly underestimated by using adiabatic results. Indeed, from our present approach we estimate the dimensionless 
spin-phonon coupling (A = ^^^) for CuGcOa one order of magnitude larger than the previously calculated one. This 
could have important consequences to interpret experimental results in this and other materials. 

A relation similar to eq. (31) holds between the transition temperature in the non-adiabatic and the adiabatic case. 
The free energy corresponding to action (20) should be expanded in powers of cq. We take again this expansion from 
Eq. (12) of [13] (now ^^'^^^2 is replaced by egg'^ F{ll>\^, ujj_) / y/a). The total free energy density is 



.J^^ _J^rp2 , 2 g'-P^(^lh^^) A _ .9^7r^f(^||,^x) \ 
■^^ NL 6vF ° 2 1^ 2r4(|)T J 



(32) 
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FIG. 1: Dependence of the gap and the SP temperature with 



The vanishing of the factor proportional to eg signals the transition temperature. We get 



and the promised relation is: 



9 9 

J^SP ^ F{uj\\,uj±_) 



Tsp 

rpad 



(33) 



(34) 



Notice that the ratio from our computation docs not depend on The value of such a ratio seems to be 
fixed by the adiabatic calculation. The existence of such universal ratio is similar to the BCS mean field theory of 
superconductivity. Moreover, the experimental value observed in CuGeOa is close to the BCS-ratio 1.76. In fact 
a BCS-like theory could be obtained from the original Heisenberg model coupled to phonons by a Jordan- Wigner 
transformation and neglecting interaction between the spinless fermions. Note however that the spinless fermions are 
strongly interacting, and this fact renormalizes the BCS ratio away from the non interacting value^'^. The precise value 
of this ratio is difficult to obtain because it is affected by logarithmic corrections induced by marginally irrelevant 
terms. Indeed, the effect of the marginal term arising from the Heisenberg interaction is expected to be small in 
CuGeOa because the second next- nearest neighbor is closer to the critical value^^, where logarithmic corrections 
vanish. However, as can be seem from action (11), logarithmic corrections are also expected from the spin-phonon 
interaction term and they do depend on the phonon frequencies. Therefore we expect that the exact value of -jA^ will 
ultimately depend on — . 



IV. DYNAMICS OF THE PHONON MODES 



The original work of CF on spin-Pcierls systems predicts the softening of a phonon mode above the transition 
temperature, whose polarization pattern corresponds to the static distortion leading to the low temperature phase 
below the transition point Tsp. The softening is due to the coupling of the phonons with the one-dimensional 
magnetism. This is in fact a quite general scenario for a structural displacive phase transition. Therefore it came as 
a surprise that the SP transition in CuGeOa shows no phonon softening^^'^°. Gros and Werner^ have reanalyzed the 
CF approach showing that, for higher enough bare phonon frequency, there is a second pole of the phonon propagator 
without softening of the active SP phonon, characterized by the presence of a central peak in the dynamical structure 
factor near the transition point. Moreover, the RPA approach underlying the CF calculation is justified if the 
SP phonon has a bare frequency considerably lower than other momentum phonons, i.e. a non-magnetic softening 
observable at high temperature relative to the SP transition temperature. This is not the situation in CuGeOa and 
application of RPA seems questionable. 

In the following we show how by extending our previous calculations to the dynamics of the phonon modes we 
can solve this puzzle. The extension of the static MF approximation on the interchain interaction to the dynamical 
correlation function is accomplished by a RPA analysis along the lines of [12.21]. In contrast to the one-dimensional 
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analysis in CF. this RPA is based on the gMasi-one-dimensionahty underlying this kind of materials, and does not 
require a preexisting mode softening. It becomes exact in the limit of an infinite number of neighboring chains. 

The dressed phonon frequencies are obtained from the poles of the retarded phonon Green's function D{k,q,ujy^ 
which can be computed from the Matsubara Green's function D{k, g, a;„), with wave vectors k in the chain direction 
and q the transverse one, by analytic continuation as D{k,q,uj) = D{k, q^icon uj). This is in fact the Fourier 
transform of the Green's function in real space and imaginary time D^~^ {x — x',r — r'), which can be found by 
functional derivation of the partition function when a current term is added: 



Z[J]= J 2?u^I?(/)^exp|-^ J dxdT[C{(fy>{x,T),u^{x,T))+gJ^{x,T)u^{x,T)]^ 
Integrating out the phonon coordinates as before we obtain the following effective action 



(35) 



2 



3, J 

The dressed phonon propagator 

(37) 



{JHx.r)}=0 



g2 SJi'{x',T')SJi{x,T) 
can be obtained from eqs. (35, 36, 37) as 

D^-''{x-x',t-t') = G^'-^{x' -x,t' -T)-g^Y2 J dxdx'drdT'G'^-' {x - x,t - t)g'^' {i' - x' ,? - t') 

hi' 

(t(: cos{V2^(j^{x,T)) :: cos{V2^(j/' {x' ,?)) :)), (38) 



where G^-J {x - x' ,t - t') ^ 6{x - x')G" (t, t') is a convenient notation for the non-interacting phonon propagator 
and T indicates Euclidean time-ordering. After Fourier transforming to [k^q,ojn) coordinates and performing the 
analytic continuation iw„ — > w we obtain the desired retarded phonon Green's function 

D{k,q,Lo) = G{q,Lu)+g^G\q,Lu)x{k,q,uj), (39) 

where G{q, cu) = -~uj'^+ui(q) ^'^'^ bare phonon Green's function and x(^j ^) stems for the Fourier transform of the 



retarded correlator of : cos(v27r(/)^) : . In order to evaluate bulk corrections arising from inter-chain coupling, we will 
calculate this correlator using a RPA approach in the interchain interaction^^. It reads 

x^^^l^-?,-) - 1 ^r^'^^t^'fl, rv (40) 

1 - g^Ginter[q,uj)x"[k,uj;l ) 

where x°(fc)'i';T') is the already known finite temperature one chain : cos{^/2tt(1>) :-correlator and Ginter is the inter- 
chain bare phonon propagator, defined in coordinate space by the subtraction 

GlnLi^ - x', T - r') = G^~i'{x ~x',T- t') - P^'G°{x ~x',t- t'). (41) 

The one chain : cos(-\/27r0) :-correlator at finite temperature has been computed elsewhere^'^'^^, being given by 

where d ~ 0.37 is a constant weakly depending of the momentum cutoff in the bosonization procedure, A = vp\k — 7r| 
is the lower edge of the two spinon continuum, and 



1 r 

Jo 



dxe^''^ {smh{x)) ^ 



1 Ffi 



4 " 2 



2x/2^r[|-i|] 



(43) 
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It will be useful to note that at low frequency one can expand 

Tx°(fc,.;T)«.xo-.Xi^+X.(^)\ (44) 

with xo 0.26 and X2 ~ 2.2. 

The Fourier transform of the interchain propagator G*"*'^'' is obtained (see Section II) as 

G-*-(g,c.) = ,\ . -T{uj), (45) 



with 



2tt 



We are interested in the temperature evolution of the SP active phonon mode, given by fc = tt and q = 0. From 
eq. (39-45) we have 

D(0, .) = ..-(l-..-T(.)x°(.,c.;r)) 

G-i + ^(l-G-iT(c.))x''(^,u;;T) 

and the phonon frequencies are given by the poles over the real axis of this expression, i.e. the roots of the following 
equation 

^ ^2 ^ _ (^2 _ iu^)T{u;))Rex°{7:, co; T) = 0. (48) 

Let us start the analysis of the consequences of the previous results by recalculating the spin-Peierls transition 
temperature. Note that, different than the previous Section, we are now coming from the high temperature phase to 
the transition point. As discussed by GW, Tgp is signaled by a macroscopic occupation of the Peierls active phonon 
mode, i.e. the transition temperature takes place when cq. (48) has a solution lj = Q. We then obtain 



^^^^ 



1 - 



(49) 



It is important to stress that eq. (49) gives the same dependency of Tsp on the microscopic parameters than eq. (33) 
in the previous Section; differences in the numerical prefactor arise from different criteria to fix the momentum cutoff. 
Note however that relation (34) is fulfilled. This results prove that our RPA procedure contains the previous MF 
approach when the static properties are considered. 

In the rest of this Section we face the question of phonon softening, discussing the calculation of the renormalized 
phonon frequencies. We use (49) to eliminate the cut-off dependent coupling from cq. (48). After straightforward 
algebra we obtain 



XO 1- , - -1 / 1- . ^ " Ur,pfiexV,c.). (50) 




This equation is the generalization of the GW results to finite values of and agrees with their corresponding eq. 

(6) when ^ ^ 1. There are two different temperature dependence regimes of the renormalized phonon frequency 

^(TrjO), as depicted in Fig. 2. For low enough ti>|| the phonon progressively softens from higher T until Tsp, where it 
drops to zero frequency. In this regime the phonon softening explains the SP transition. For higher W|| , the phonon 
frequency remains finite up to the SP temperature and no soft phonon behavior is obtained; in fact, there can be 
partial softening or even hardening of the renormalized frequency close to Tsp. However, an additional solution lo = 
in eq. (50) signals the appearance of a central peak in the dynamical structure factor, leading to the SP transition^. 

The limiting value of uj\\/Tsp for separation of regimes depends on lo^. To find this transition analytically it is 
enough to compare concavities of both sides of eq. (50) at T = Tsp, as seen in Fig. 2. Expanding both sides up 
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FIG. 2: Evolution of the SP relevant phonon mode. Dashed lines, from blue to red, represent the rhs of eq. (50) when 
temperature is lowered until Tsp. Solid gray lines represent the Ihs for lj^^/Tsp — 2,3,4 and uij_/Tsp = 1. The intersections 
illustrate the evolution of the renormalized phonon frequency. For ui^^/Tsp — 2 the phonon softens completely, for ui^^/Tsp = 3 
it softens partially, and for ijj\\/Tsp = 4 it hardens close to Tsp- 



to second order in the soft phonon regime is reached when, at T = Tsp, the coefficient of ijP in the r.h.s. is 

smaller than the one in the l.h.s. Under this condition the bare phonon frequency (that obtained at hight T) evolves 
continuously to zero when the temperature is lowered until the transition temperature. The following inequality 
should then be fulfilled to have the soft phonon regime: 



-11 iUi^^TT^U^^ ^^^^ 



Tsp V X2 V 2 1 + {uj±/u:\\Y 



In the other case, the bare phonon frequency evolves to some finite value at the SP transition. In Fig. 3 we show 
the phase diagram of the spin-phonon system on the uj±^ — lo\\ plane. Note that in a wide range of high enough 

the separatrix is given by ^ 2.2 as in GW. Our result extends this condition for finite values of where 
realistic spin-Peierls materials as tJuGeOa live. This also explains why the RPA approach for the spin-phonon coupling 
underlying the GW calculation compares well with experiments on CuGeOa. 

In Fig. 4 we show the temperature dependence of the renormalized phonon frequency for ~ 1 where the 
phase diagram of Fig. 3 is not completely flat. The behavior shifts from a soft to a non soft phonon regime at some 
intermediate value of ;f^- In this last regime the mechanism of the SP transition has been associated with the 
emergence of a central peak in the spectral signal measured in neutron scatcring experiments^^. For high enough 

the phonon hardens from a minimum value at some intermediate temperature. kX, T = Tgp the frequency of 
the phonon is the bare one as in the high temperature limit. We will use this fact to fix uj^ for CuGeOa. Note that 
the T-dependence of the phonon is similar to the one obtained by GW but the frequency and temperature scale are 
different. 

From the previous results we can follow the T-dependence of the Peierls active phonons in CuGeOa. In Fig. 4 (6), 
we show our prediction for the higher energy mode T^, which is the most anomalous and strongly coupled to the 
magnetism"*. We take as a bare phonon frequency lo\\ ~ 6.8THz i.e. its low temperature value. For uj± we use the 
relationship obtained in the previous Section. Note that a value of Tsp = 14.1ii'(0.294TiJz) has been used to give 
the frequencies in Thz and the temperatures in K. Our results compare well with Fig. 3 of [4]. In fact we predict a 
stronger hardening than the one experimentally observed (7% against 4.5% seen in the experiment). This could be 
due to the simplified model we use, which does not take into account the fact that two phonon modes are necessary 
to describe the SP transition in CuGcOa"'. 

Besides, we can make a prediction for the phonon frequency behavior at higher temperature than the one already 
measured. It should increase as shown in Fig. 4, approaching the bare frequency at very high T. We note that at 
room temperature the phonon frequency does not correspond to the bare one as expected. It is almost the minimum 
of the predicted phonon frequency. The bare frequency is only obtained at T ^ Tsp and at very high temperature 
(probably not experimentally accessible). 

To summarize the results of the present Section, we have generalized the MF approach of the previous Section to 
the dynamical correlation functions in the high temperature phase. Our approach relies on a RPA on the effective 
interchain coupling generated by the spin-phonon interaction. In contrast to the previous theoretical calculations. 
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FIG. 3; The line separating the zone in the t^^, plane, where the spin-Peierls transition takes place by softening of the 
phonon from the one where the Peierls active phonon does not soften completely. 



03, =3.94 THz m,,=6.8 THz 




FIG. 4: (a)The temperature dependence of the renormalized phonon frequency for different values of the bare with 

(b) The hardening of the phonon for the parameters describing the higher energy active SP mode in CuGeOs. The bare values 

used for clI|| and cjx are shown on top of the Figure. 



the present one does not presuppose the existence of a soft phonon mode in the non magnetic regime. Therefore, it 
can been apphed to CuGeOa. Moreover, the quaUtative behavior of the Peierls active phonon is similar to the one 
previously obtained. There is a regime corresponding to low frequency phonon where the structural phase transition 
takes place by progressive softening of a phonon from the high temperature phase. For higher bare phonon frequency 
no soft phonon is observed and it can even harden from the room temperature value to the transition temperature 
one. This is precisely the behavior seen in neutron scattering measurements of the phonon spectra in CuGeOa. 



V. CONCLUSIONS 



We have developed a route for the spin-Peierls transition which goes beyond the usual adiabatic treatment of 
the phonon field. We emphasized the essential character that plays the dispersion of the phonon in the transversal 
direction to the magnetic chains. Moreover, the results of the present paper show that, when the adiabatic hypothesis 
for the phonon coordinates is relaxed, the one-chain model does not represent a good starting point to describe a 
real system. Furthermore the effective interchain interaction generated by the phonons is an essential ingredient for 
the opening of the gap and the existence of a finite temperature phase transition. The precise determination of the 
SP transition temperature from the microscopic parameters strongly depends on the width of the transversal phonon 
dispersion. 

Our approach for the dynamics of the phonon mode in the high temperature phase justifies and generalizes previous 
calculations of the phonon spectra based on the Cross and Fisher early work. Applied to CuGeOa, we show that 
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the phonon does not soften until the phase transition. Moreover the overall temperature dependence of the phonon 
frequency is consistent with experimental determinations. 

As the elastic interaction between the chains play a central rol in non-adiabatic spin-Peierls systems, we anticipate 
that a change in the relative positions of the magnetic ions in different chains will change the phase diagram and 
the mechanism of the phase transition. In this sense it is highly interesting to study the recently discovered quasi- 
one-dimensional magneto-elastic system TiOCl with an intermediate incommensurate phase where the Ti ions are 
not aligned in the direction perpendicular to the magnetic chain'^. We will study in a forthcoming work this system 
within the formalism developed in the present paper. 
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